Universal spin-Hall conductance fluctuations in two dimensions 
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We report a theoretical investigation on spin-Hall conductance fluctuation of disordered four 
terminal devices in the presence of Rashba or/and Dresselhaus spin-orbital interactions in two 
dimensions. As a function of disorder, the spin-Hall conductance GsH shows ballistic, diffusive and 
insulating transport regimes. For given spin-orbit interactions, a universal spin-Hall conductance 
fluctuation (USCF) is found in the diffusive regime. The value of the USCF depends on the spin- 
orbit coupling tso, but is independent of other system parameters. It is also independent of whether 
Rashba or Dresselhaus or both spin-orbital interactions are present. When tao is comparable to the 
hopping energy t, the USCF is a universal number ~ 0.18e/47r. The distribution of Gsh crosses over 
from a Gaussian distribution in the metallic regime to a non-Gaussian distribution in the insulating 
regime as the disorder strength is increased. 

PACS numbers: 71.70.Ej, 72.15.Rn, 72.25.-b 



The notion of dissipation- less spin-current^ has at- 
tracted considerable interests recently. In its simplest 
form, a spin-current is about the flow of spin-up elec- 
trons in one direction, say ~\-x, accompanied by the flow 
of equal number of spin-down electrons in the opposite 
direction, —x. The total charge current in the a;-direction 
is therefore zero, /g = e(/| -I- J|) = 0; and the total spin- 
current is finite: Ig = h/2{I^ — I^) ^ 0. For a pure semi- 
conductor system with spin-orbital (SO) interactions, it 
has been shown^ that an electric field in the z-direction 
can induce the flow of a spin-current in the x-direction 
perpendicular to the electric field: such a spin-current is 
dissipation-less because the external electric field does no 
work to the electrons flowing inside the spin-current. If 
the semiconductor sample has a finite a;-extent, the flow 
of spin-current should cause a spin accumulation at the 
edges of the sample, resulting to a situation that spin-up 
electrons accumulate at one edge while spin-down elec- 
trons at the opposite edge. Hence a spin-Hall effecfeS*^ 
is produced where chemical potentials for the two spin 
channels become different at the two edges of the sam- 
ple. This interesting phenomenon has been subjected to 
extensive studies and there are several experiments re- 
porting spin accumulation which may have provided ev- 
idence of this effecti. It has been shown that for a pure 
two dimensional (2d) sample without any impurities, the 
Rashba SO interaction generates a spin-Hall conductiv- 
ity having universal valued of e/87r. It has also been 
shown that any presence of weak disorder destroys this 
spin-Hall effect in the large sample limit^. On the other 
hand, numerical studies have provided evidence that for 
mesoscopic samples, spin-Hall conductance can survive 
weak disorde r ''i^i^ . 

One of the most striking quantum transport features 
in mesoscopic regime is the universal charge conduc- 
tance fluctuation (UCF)i24iiiiS: quantum interference 
gives rise to the sample-to-sample fluctuation of charge 
conductance of order /h, independent of the details of 



the disorder, Fermi energy, and the sample size as long as 
transport is in the coherent diffusive regime characterized 
by the relation between relevant length scales, I < L < 
Here L is the linear sample size, I the elastic mean free 
path and ^ the phase coherence length. If time-reversal 
symmetry is broken, UCF is suppressed by a factor of 
two. 



A very important and interesting issue therefore arises: 
what are the properties of the fluctuations of spin-Hall 
conductance in disordered samples? Is there a transport 
regime where spin-Hall conductance fluctuation is uni- 
versall The answers to these questions are non-trivial 
because the flow of dissipation-less spin-current is quali- 
tatively different from the flow of charge current driven 
by an external electric fleld. It is the purpose of this 
paper to report our investigations of these issues. For 
a disordered four-terminal sample with a given Rashba 
SO interaction strength tgo, and/or Dresselhaus interac- 
tion strength tso2, our results suggest that there is indeed 
a universal spin-Hall conductance fluctuation (USCF) 
whose root mean square amplitude m g — 0.18(e/47r), in- 
dependent of other system details (thus universal). The 
fluctuation is however a function of the SO interaction 
strength and found to be well fitted by a functional form 
of rms(Gsi/) = gtanh(|tso — iso2|/0.17). Finally, the dis- 
tribution of spin-Hall conductance obeys a Gaussian dis- 
tribution in the metallic regime and deviates from it in 
the insulating regime. 



To investigate USCF, we consider a four terminal de- 
vice in two dimensions schematically shown in the left 
inset of Fig^ We will first discuss the results in the 
presence of only Rashba interaction. In the presence of 
Rashba interaction (a^oZ ■ (c x k)), the Hamiltonian of 
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FIG. 1: Spin-Hall conductance fluctuation versus energy 
for disordered samples. Triangles, squares, and circles are 
for W = 1,2,3 respectively. Left inset: schematic plot of 
the four-terminal mesoscopic sample where Rashba interac- 
tion exists in the center scattering region and the leads 1, 3. 
The width of the square sample is L. A small voltage bias 
is across leads 1,3, and spin-Hall conductance is measured 
through leads 2,4. Right inset: the ensemble averaged spin- 
Hall conductance GsH versus electron energy for SO interac- 
tion strength tso ~ 0.3. Solid line: pure sample with W = 
and other symbols are the same as the main panel. In all 
figures the spin conductance and its fluctuation are measured 
in e/47r. 



this device is given hy^: 

H = -t J2 ('=L"=J'^ + ^■^■) + J2 '^^4a(^^<^ 

<ij>(T i<7 

i 

-i(4^Q + y,X c\^C^ + y^^) + h.c] (1) 

where cj^. is the creation operator for an electron with 
spin cr on site i, x and y are unit vectors along x and 
y directions. Here t = /2ma? is the hopping en- 
ergy and tso = OLsol^a is the effective spin-orbit cou- 
pling. The on-site energy is given by = 4t. In addi- 
tion, static Anderson-type disorder is added to with a 
uniform distribution in the interval [—H^/2, 14^/2] where 
W characterizes the strength of the disorder. We con- 
sider the situation where Rashba interaction is present 
everywhere except in leads 2 and 4 (see the left inset of 
FigD in order to measure the conserved spin-current^. 
We apply external bias voltages at the four terminals as 
(Vi,i = 1---4) = (w/2, 0, — w/2, 0): such a setup gen- 
erates a spin-current flowing from lead 2 to 4, i.e. a 
spin-Hall effect measured from these two leads* 

The spin Hall conductance GsH is calculated from the 



Landauer-Biittiker formula™ 

GsH = (e/87r)[(r2Ta - ^24,1) - {T^^.z ~ T^x.z)\ (2) 

where the transmission coefficient is given by T^a^x = 
Tr(r2o-GTiG'°), here C'" are the retarded and advanced 
Green's functions of central disordered region of the de- 
vice which we evaluate numerically. The quantities Vi„ 
are the line width functions describing coupling of the 
leads to the scattering region, and are obtained by calcu- 
lating self-energies due to the semi-infinite leads using a 
transfer matrices methodic. The spin-Hall conductance 

fluctuation is defined as rms(Gs/f) = ^ {G^fj) — (Gsh)^, 
where (■ ■ •) denotes averaging over an ensemble of sam- 
ples with different configurations of the same disorder 
strength W. Note that in the presence of disorder, 
although one could use another the definition GsH — 
(e/47r)(T2f _i — T2i^i) to calculate and obtain the same 
average spin-Hall conductance as that of Gsh, the spin- 
Hall fiuctuation can only be obtained correctly using the 
definition of Eq. (|2Jl . We perform our calculations on LxL 
square samples with four leads described above. Sample 
sizes of L = 40 up to 100 are examinediS:. To fix units, 
throughout this paper we measure the energy E, disorder 
strength W , and spin-orbit coupling tso in terms of the 
hopping energy t. 

Figdplots spin Hall conductance fluctuation vs Fermi 
energy at a fixed tso = 0.3 and sample size L = 40 for 
several disorder strengths W = 1,2,3. We have also 
shown the averaged spin Hall conductance in the right 
inset. Due to the electron/hole symmetry, only data for 
energy range [0, 4] are shown. Over 10 000 disordered 
samples are averaged. For pure sample without disor- 
der, as the Fermi energy is increased the number of sub- 
bands increase. As a result, the spin-Hall conductance 
GsH shows small oscillations. When disorder is increased 
from zero, GsH decreases as expected and eventually, the 
small oscillation due to the subbands vanishes. Most im- 
portantly, Fig^ shows substantial sample to sample fiuc- 
tuations of GsH, measured by Tms{GsH), of the order 
6^ where (5 is a number between 0.1 and 0.2. Such an 
amplitude of fluctuation is comparable to the spin-Hall 
conductance itself. 

In Fig. 2a, b, we plot the (Gsh) and rms(Gs/f) as a 
function of disorder strength W at fixed Fermi energy 
E = I for a number of different spin-orbit couplings 
from tso — 0.2 up to 0.7. Several observations are 
in order. First, the spin-Hall conductance (Fig[2jL) de- 
creases smoothly as the disorder strength is increased: 
the transport characteristics goes from quasi-ballistic at 
small W to the diffusive regime at larger W. In the dif- 
fusive regime, the spin-Hall conductance decreases expo- 
nentially with the disorder strength between W = [1, 5]. 
Finally it goes to the insulating regime for even larger 
W where GsH vanishes. Second, the numerical data 
show that the onsets of insulating regime Wc are dif- 
ferent for different spin-orbit couplings tso- The larger 
the spin-orbit coupling, the larger Wc- This finding is 
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FIG. 2: (a) Ensemble averaged Gsh versus disorder strength 
W for tso ~ 0.2 (cross), 0.3 (filled triangle), 0.4 (open circle), 
0.5 (star), 0.6 (rhombus), 0.7 (filled square). The average is 
over 20 000 samples with L = 40. Inset: size dependence of 
spin-Hall conductance fluctuation with tso ~ 0.3. Different 
symbols are for W = 1 (stars), 2 (rectangles), 3 (circles) and 
4 (triangles). The ensemble average is over 20 000 samples 
for different size L. (b) The corresponding ensemble averaged 
spin-Hall conductance fluctuation versus W , the symbols are 
for the same tso values as in (a). 



consistent with that of Ref 8 which suggested that the 
localization length depends on tso and belongs to two- 
parameter scaling. Third, from the spin-Hall conduc- 
tance fluctuation shown in FiglSb, we observe that for 
small disorder W < 1, the fluctuations for different tgo 
have very similar values and the curves collapse. For 
larger disorder, the fluctuations develop a plateau struc- 
ture so that tt[is{Gsh) becomes independent of the dis- 
order parameter W for each given tgo- In this sense, the 
fluctuation rms(Gs_f/) becomes "universal" and the spin- 
Hall transport enters the regime with universal spin-Hall 
conductance fluctuations. Importantly, both the width of 
the plateau and the value of USCF depend on tso- The 
larger the tso^ the wider the fluctuation plateau which 
characterizes the diffusive regime. 

Now we examine the dependence of spin-Hall conduc- 
tance fluctuation rms(G's_f/) on system size L in the inset 
of FigEt for tso = 0.3. With weak disorder W ^ 1 
(stars) the fluctuation increases with sample size indicat- 
ing that spin-Hall conductance is not yet in the USCF 
regime because transport is quasi-ballistic. In the dif- 
fusive regime, W = 2,3,4, the fluctuations saturate at 
i5e/47r where S ~ 0.2. The independence of system 
size by the fluctuation Tms{GsH) provides a strong evi- 
dence of USCF. Namely, as long as transport is in the 
diffusive regime, the fluctuation of the spin-Hall con- 
ductance is dominated by quantum interference giving 
rise to a universal amplitude. Since the value of USCF 
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FIG. 3: (a). USCF values versus Rashba spin-orbit coupling 
tso at = 1. Inset: rmsGsH versus Fermi energy for W — 3 
in the presence of both Rashba and Dresselhuss SO coupling, 
tso = 0.5 and tso2 ~ 0.2. (b)/(c). rmsGsH versus Desselhaus 
SO coupling tso2 at B = 2, W = 3, and tso = 0.5/tso = 0.1. 
In the cases of inset of (a), (b) and (c), 40000 samples are 
collected. 



(which is obtained from each curve in the lower panel 
of Fig. 2) depends on the spin-orbit couphng tso as seen 
from Fig. 2, we have obtained a collection of the USCF 
for different tso which is shown in Fig. 3a. Interestingly, 
the USCF can be well fltted (solid line) by a function 
rms(GsH) ^ 5tanh(tso/0.17) where g = OASe/An. This 
can be understood as follows. When tgo = 0, there is no 
spin-Hall current and hence no fluctuations. There is a 
crossover regime before the fluctuation saturates to the 
USCF plateau. 

Fig0Ji-d plot the distribution function of spin Hall 
conductance, P{Gsh), for several different disorder val- 
ues W. For each W, data is accumulated by calcu- 
lating 20 000 reahzations of disorder. P{Gsh) appears 
to clearly obey a Gaussian distribution in the metallic 
regime up to « 5. For larger W between [5, 10], trans- 
port is in the insulating regime, the symmetric distribu- 
tion exhibits non-Gaussian behavior (Fig^). At even 
larger disorder W — 12 shown in Fig^, the distribution 
becomes non- Gaussian and asymmetric. The deviation 
from Gaussian distribution can be characterized by the 
moments of spin-Hall conductance. We have calculated 
the skewness 71 and kurtosis 72 whose definitions arei^, 

71 = ^J.3/^A^'^ and 72 = Ai4//i|-3 where /i„ = ((a;- (a;))") 
(n — 2,3,4) denote the central moments. The skew- 
ness describes the degree of asymmetry of a distribution 
around its mean while the kurtosis measures the relative 
peakedness of a distribution. The results are plotted in 
Fig2f, showing that in the metallic regime Vl^ < 5, both 
skewness and kurtosis are essentially zero while they be- 
come non-zero for larger W, consistent with the distri- 
butions. Hence the skewness and kurtosis can be used to 
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FIG. 4: (a-d) The distribution of the spin-HaU conductance 
for different disorder strengths at a fixed energy E = 1 and 
tso = 0.3. Data are collected for 20 000 samples for each W 
with L = 40. (e) The skewness 71 and kurtosis 72 versus dis- 
order strength W for the same ensemble, (f) rmsGsH versus 
Fermi energy for W = 3 and tso = 0.6, 0.7. 



identify the diffusive regime. Importantly, these quanti- 
ties can be measured experimentllyi^. Finally, we have 
checked that the above features of the spin-Hall conduc- 
tance fluctuation are generic and valid for other values of 
E and tso- For instance. Fig. 4f shows the rmsGsH versus 
Fermi energy when W = 3 and tso — 0.6, 0.7. We see 
that between E = I and E = 3, rmsGsH is around the 
universal value 0.18. 

So far we have focused on spin-Hall conductance fluc- 
tuations with the Rashba interaction. To further demon- 
strate the universal behavior, we have also analyzed the 
case of Dresselhaus spin-orbital interaction by adding a 
term l3so{o'xkx — (^yky) in Eq.(l). Using unitary trans- 
formations, it is easy to prove that for the spin-Hall 
current along z-direction we have /f^(aso = Q,Pso) ~ 
IsH{oi'>o,l3so = 0) and I^Hic^so = Pso) = 0. When 
there is no Rashba spin-orbital interaction (ago — 0) and 
only Dresselhaus term exists, we have obtained exactly 
the same USCF value and behavior as that of Rashba 
term alone. When both Rashba and Dresselhaus terms 
are present, it is not obvious that the USCF persists. 



This is because these two terms have different symme- 
try, the Rashba coupling arises from the structure inver- 
sion asymmetry with SU(2) symmetry while the Dres- 
selhaus coupling arises from the bulk inversion asym- 
metry with SU(1,1) symmetrjii^. From Fig. 3b and 3c 
plot we see that USCF for the latter situation (both 
Qfso and Pso 7^ 0), the results are similar to the case 
of pure Rashba or pure Dresselhaus interaction. Because 



Pso) 



0, the USCF curves have a dip to 



zero when a^o — Pso- Defining tso2 = Pso/'^a, Figs. 3b, 
3c show that USCF is reached when \ts02 — tso\ ^ 0.4. 
Finally, the inset of Fig. 3a shows that the USCF is inde- 
pendent of Fermi energy when both SO interactions are 
present. 

In summary, for spin-Hall effect generated by the 
Rashba and Dresselhaus interactions in mesoscopic sam- 
ples, our results strongly suggest the existence of a uni- 
versal spin-Hall conductance fluctuation due to impu- 
rity scattering in the quantum coherent regime. In this 
regime, the USCF is characterized by sample-to-sample 
fluctuations of GsH for a given SO interaction (Rashba 
or/and Dresselhaus) strength tso^ measured by quantity 
rms(Gsij) with a universal amplitude g = 5-^ where 
5 ~ 0.18, which is independent of system size, impu- 
rity scattering strength, and Fermi energy. Importantly, 
this fluctuation amplitude is of the same order as the 
spin-Hall conductance itself. Comparing with the fa- 
miliar UCF in charge conductance of disordered meso- 
scopic samples, USCF originates from a similar physics 
of quantum interference effect which leads to significant 
sample-to-sample fluctuations in spin-Hall conductance. 
A main difference is that spin-Hall effect is due to SO 
interactions in the sample, thereby the USCF also de- 
pends on the SO parameter tgo (or ^502)) and our nu- 
merical results can be well fitted by a functional form of 
rms(Gs//) = gtauh{\tso ~ tso2|/0.17). 
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